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ABSTRACT 


The  laminar,  single-phase  heat  transfer  and  friction  characteristics  of  a porcupine-like  surface  (integral- 
spine-fin)  within  an  annulus  are  presented.  The  heat-transfer  coefficient  was  determined  using  a modified 
version  of  the  Wilson  Plot  method  on  a 3 m test  section.  Three  fluids  were  investigated:  (1)  water,  (2) 
34%  ethylene  glycol/  water  mixture,  and  (3)  40%  ethylene  glycol/water  mixture.  These  fluids  produced 
a significant  variation  in  the  Prandtl  number  so  that  its  exponential  dependence  could  be  determined.  The 
annulus  Reynolds  numbers  were  varied  from  100  to  1400  to  obtain  the  Reynolds  number  exponent.  An 
empirical  correlation  for  the  Nusselt  number  that  accounts  for  the  development  of  the  thermal  boundary 
layer  is  presented.  An  empirical  correlation  for  the  fanning  friction  factor  is  also  provided.  It  is  shown 
that  the  spines  enhance  the  heat  transfer  through  additional  surface  area  and  fluid  mixing. 

Keywords:  annulus,  enhanced  heat  transfer,  spine-fin,  friction,  water,  ethylene  glycol 
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NOMENCLATURE 
English  Symbols 

A surface  area  utilized  in  heat  transfer  (m2) 

A„  cross  sectional  flow  area  (m2) 

A0  total  surface  area:  Ar  + Af  (nr) 
a tube-side  Reynolds  number  exponent 

b annulus-side  Reynolds  number  exponent 

C empirical  constant  in  Nusselt  number  correlation 
cp  specific  heat  (J/kg»K) 

D0  tube  diameter  over  the  fins  (m) 

Dr  tube  diameter  at  the  base  of  the  fins  (m) 

Dp  pre-finning  pipe  outer  diameter  (m) 

Dh  hydraulic  diameter  of  the  annulus  side:  4Ac/pw  (m) 
Di  inner  diameter  of  the  pipe  (m) 

f fanning  friction  factor 

Gr  Grashof  number:  Dh3p12/S(TW  - Tm)//* 

h heat-transfer  coefficient  (W/m2  • K) 

k thermal  conductivity  (W/m  • K) 

L axial  length  of  tube  (m) 

m mass  flow  rate  (kg/s) 

Nf  number  of  fins  per  axial  length  (m'1) 

Nu  Nusselt  number 

n annulus  Prandtl  number  exponent 

P pressure  (Pa) 

Pr  Prandtl  number 

p exponent  on  thermal  developing  correction 

pt  axial  fin  pitch  (m) 

pw  wetted  perimeter:  7rD0  + A0/L  (m) 
q heat  load  (W) 

q"  heat  flux  (W/m2) 

Rea  annulus  Reynolds  number:  4rha/(/4a(7rD0  + A0/L)) 
Ret  smooth  tube  Reynolds  number:  4mt/(/zt7rDi) 
s estimate  of  standard  deviation 

T fluid  temperature  (K) 

Tw  temperature  of  spine-fin  tube  wall  at  root  of  fin  (K) 

t thickness  of  fin  (m) 

UA  overall  conductance  (W/K) 

X lumped  parameter  on  Wilson  plot  ordinate 
x thermal  entry  length  (m) 

Y lumped  parameter  on  Wilson  plot  abscissa 


Greek  symbols 

/3  coefficient  of  volumetric  thermal  expansion  (K'1) 

AP  pressure  drop  (Pa) 

ATlm  log-mean  temperature  difference  (K) 

Ax  distance  between  pressure  taps  (m) 

77  surface  efficiency 


?7f  fin  efficiency 

e (Tw  - Ti)/(TW  - Tm) 

p dynamic  viscosity  (kg/m  • s) 

p density  (kg/m3) 


Subscripts 

a annulus-side  fluid 
cu  copper 

f fin 

i inlet 

1 liquid 

m mean,  bulk 

r root 

t tube-side  fluid 


INTRODUCTION 


This  paper  presents  empirical  correlations  for  the  single-phase  Nusselt  number  and  fanning  friction  factor 
for  an  integral-spine-fin  within  an  annulus.  A modified  Wilson  Plot  was  used  to  derive  the  Nusselt 
numbers  from  overall  conductance  measurements.  The  friction  factors  were  calculated  from  pressure 
drop  measurements.  The  test  fluids  --  water,  34%  ethylene  glycol/water  mixture,  and  a 40%  ethylene 
glycol/water  mixture  --  were  chosen  to  provide  a significant  variation  in  the  Prandtl  Number  (4  to  22). 
In  this  way,  the  exponential  dependence  of  the  Prandtl  number  on  the  heat  transfer  was  experimentally 
determined.  For  an  analogous  reason,  the  annulus  Reynolds  numbers  were  varied  from  100  to  1400. 
Average  heat  fluxes  to  the  test  section  ranged  from  13.5  to  54  kW/m2. 

Very  little  friction  and  heat  transfer  data  exist  for  the  integral-spine-fin.  Most  studies  have  been  done 
for  tubes  in  cross  flow  (Holtzapple  and  Carranza  (1990)).  The  work  of  Hobson  and  Weber  (1954)  is, 
to  the  authors’  knowledge,  the  only  reported  work  to  use  the  integral-spine-fin  within  an  annulus. 
However,  the  fin-tips  in  their  study  did  not  touch  the  outer  shell  which  allowed  the  fluid  to  bypass  the 
fins.  Also,  the  flow  was  turbulent.  In  contrast,  all  of  the  annulus-side  Reynolds  numbers  in  the  present 
study  are  laminar  and  the  fins  touch  the  outer  shell. 

Both  radial  and  axial  cross  sectional  schematics  of  the  integral-spine-fin  tube  are  given  in  Figure  1.  The 
tube  has  approximately  8189  fins  per  meter  (Nf)  and  an  axial  fin  pitch  (pt)  of  3.175  mm.  The  diameter 
of  the  tube  over  the  fins  (D0)  is  26.8  mm.  The  tube  diameter  at  the  base  of  the  fins  (Dr)  is  11.72  mm. 
The  root  surface  area  per  unit  length  --  A^L  = 32.76  mm  — was  calculated  from:  7rDr(l  - tjVpt),  where 
an  average  tr  of  0.35  mm  was  used.  The  hydraulic  diameter  of  the  annulus  side  with  the  fin-tips  snugly 
touching  the  outer  tube  wall  is  4.59  mm.  The  inner  diameter  of  the  pipe  (Dj)  is  9.53  mm.  The  fins  are 
not  quite  90°  with  respect  to  the  tube  axis.  However,  several  tests  were  conducted  which  showed  that 
the  flow  direction  with  respect  to  the  fin  angle  does  not  have  a measurable  effect  on  the  heat  transfer  for 
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the  present  flow  conditions. 


The  integral -spine-fins  are  formed  by  lifting  the  metal  from  the  outer  wall  of  a smooth  copper  pipe.  In 
this  study,  the  pre-finning  pipe  outer  diameter  (Dp)  was  13.7  mm.  For  a given  outer  tube,  finning  of  the 
copper  pipe  results  in  6.48  times  the  surface  area  and  the  same  cross  sectional  (Ac)  flow  area  as  the  pre- 
finning  geometry.  As  shown  in  Figure  2,  the  spines  of  our  tube  resembled  a copper  tadpole  with  its 
twisted  tail  integrally  attached  to  the  root  surface  of  the  tube.  The  spines  of  different  tubes  do  not 
necessarily  resemble  tadpoles.  The  spine  fin  cross  section  may  vary  depending  on  the  number  of  spines 
per  turn,  the  base  material,  and  the  condition  of  the  tool  used  to  produce  them  (Vandersip,  1993).  The 
tail  cross  section  of  the  spine  in  this  study  was  nearly  triangular.  The  cross  section  of  the  spine  head  or 
fin-tip  was  trapezoidal.  The  surface  area  of  one  fin  was  estimated  from  the  measurements  shown  in 
Figure  2 to  be  approximately  3 x 10'5  m.  The  total  outside  surface  area  per  unit  length  --  A0/L  = A^L 
4-  NfAf  — was  estimated  to  be  0.279  m.  The  volume  of  metal  in  all  of  the  spines  calculated  using  Figure 
2 and  that  calculated  from  the  reduction  in  the  pipe  wall  after  finning  agreed  within  1%. 


If  measurements  of  the  spine’s  dimensions  are  unavailable,  the  total  surface  area  (root  and  fin  surface 
area)  of  other  spine  fin  geometries  may  be  estimated  from: 


^2  = nD 
L 1 


7t  Dztz 
Pt 


0 . 4 mm 


- Dl) 


0 . 4mm  Nf  (D0  - Dr) 


(D.  - Dr)  Nf  <!> 


The  above  equation  approximates  the  fin  cross  section  as  rectangular  and  assumes  that  one  of  the  sides 
is  0.4  mm  wide.  The  width  of  the  other  spine  side  was  obtained  by  setting  the  volume  of  metal  lifted 
from  the  tube  to  form  one  fin  equal  to  the  volume  of  a rectangular  spine. 
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APPARATUS  AND  FLOW  CONDITION 


Figure  3 shows  a schematic  of  the  test  rig  which  was  used  to  simultaneously  obtain  heat  transfer  and 
friction  data.  A cross  section  of  the  straight  3 m test  section  is  shown  in  Fig  3.  The  test  fluid  was 
pumped  through  the  annulus.  Tap  water  was  pumped  through  the  smooth  inner  tube  at  a steady  flow  rate 
of  approximately  0.18  kg/s  to  ensure  large  turbulent  Reynolds  numbers.  The  mass  flow  rates  were 
measured  with  a coriolis  flow  meter  having  an  estimated  uncertainty  of  0.001  kg/s  for  a 99.7% 
confidence  interval.  (All  uncertainties  are  for  a 99.7%  confidence  interval  unless  otherwise  stated.)  The 
temperature  difference  (AT)  for  each  fluid  stream  was  measured  with  a ten  element  thermopile  to  within 
0.01  K.  Mixers  were  used  before  the  thermopiles  to  ensure  the  measurement  of  bulk  temperatures.  The 
condenser  and  evaporator  of  a refrigeration  system  (shown  as  heat  exchangers  in  Figure  3)  were  used  to 
supply  and  remove,  respectively,  the  heat  of  the  two  fluid  streams.  The  condenser  and  evaporator 
streams  were  switched  to  achieve  both  annulus-side  heating  and  cooling  conditions.  Two  differential 
pressure  transducers,  accurate  to  within  0.26  kPa,  were  used  to  measure  the  pressure  drop  of  each 
stream.  Table  1 shows  the  estimated  percent  uncertainty  for  various  measurements.  The  uncertainty 
estimates  given  in  Table  1 are  the  root-sum-square  of  the  component  uncertainties. 

As  shown  in  Figure  3,  the  pressure  taps  were  located  far  enough  from  the  tube  bends  to  ensure  that  the 
pressure  drop  measurements  were  for  fully  developed  hydraulic  flow.  Mattingly  and  Yeh  (1991)  have 
presented  a method  for  estimating  the  decay  of  turbulent  swirl  induced  by  a 90°  bend.  The  swirl  effect 
increases  for  increasing  Reynolds  Number.  Twenty-four  tube  diameters  are  required  for  swirl  to  dissipate 
to  1%  of  its  greatest  value  for  Ret  = 40,000  (our  largest  Ret).  On  the  smooth  tube  side,  we 
conservatively  used  32  diameters  so  that  the  pressure  drop  measurements  would  not  be  affected  by  the 
90°  bend.  On  the  annulus  side,  the  effect  of  hydraulically  developing  flow  and  90°  bends  on  the  pressure 
drop  measurements  should  be  minimal  since  the  taps  were  located  66Dh  from  the  bends,  and  mixing 
induced  by  the  fins  should  quickly  establish  the  flow. 
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Ideally,  the  tube-side  heat-transfer  coefficient  should  be  nearly  constant  along  the  tube  length  for  accurate 
Wilson  plot  results.  A nearly  constant  heat-transfer  coefficient  can  be  achieved  if  the  flow  is  thermally 
developed.  Deissler  (1955)  shows  that  the  thermal  entry  length  is  approximately  3 tube  diameters  for 
turbulent  flow  inside  a smooth  circular  tube.  If  the  solution  given  by  Deissler  (1955)  is  integrated  over 
the  entrance  region,  the  average  tube-side  heat-transfer  coefficient  for  the  present  test  section  would  be 
less  than  1%  greater  than  the  constant,  thermally  developed  value.  Consequently,  the  tube-side  entrance 
region  for  the  present  test  rig  should  have  a negligible  effect  on  the  accuracy  of  the  Wilson  Plot  method. 

Conversely,  the  thermal  entry  length  (x)  on  the  annulus  side  is  significant  for  most  of  the  171  heat 
transfer  data  points.  It  was  estimated  from  an  expression  — x = 0.05DhRePr  --  given  by  Kays  and 
Crawford  (1980)  to  be  from  16%  to  66%  (38%,  on  average)  of  the  3m  test  section.  Consequently,  this 
paper  empirically  correlates  thermally  developing  flow  data.  Considering  that  laminar  flow  exists  for 
many  liquid  heat  exchangers  in  ground  source  heat  pumps  and  other  equipment,  a predictive  correlation 
should  include  a correction  for  thermally  developing  flow.  Although  the  local  heat  transfer  varies  along 
the  entry  length,  the  accuracy  of  the  Wilson  plot  method  is  not  compromised  since  the  variation  is 
accounted  for  by  using  an  appropriate  expression  in  the  correlation  for  the  average  Nusselt  number. 


FANNING  FRICTION  FACTOR 


Mass  flow  rate  (m)  measurements  and  differential  pressure  (Ap)  measurements  across  a 2.47  m length 
of  the  annulus  (Axa)  and  a 3.46  m length  of  the  inner  tube  (Axt)  were  used  to  calculate  the  fanning 
friction  factors  for  the  respective  flows.  The  fanning  friction  factor  for  the  smooth  tube  side  was 
calculated  from: 


-Pi  APc7t2r>l 

32  fnt2 Axt 
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The  uncertainty  of  ft  is  estimated  to  be  1.5%  of  ft. 


The  annulus-side  fanning  friction  factor  was  calculated  from: 


f 


a 


2P  i.VAPa 

Axa  fna2(izD0  + 


(3) 


The  uncertainty  of  fa  is  estimated  to  be  0.003  kg/s.  That  is,  the  uncertainty  is  10%  and  2%  of  fa  for  the 
lowest  flow  rate  (0.03  kg/s)  and  the  highest  flow  rate  (0.15  kg/s),  respectively. 

Two-hundred-and-eighty-one  data  points  were  used  to  fit  the  measured  friction  factors  to  the  measured 
Reynolds  number.  The  smooth-tube-side  fanning  friction  factor  correlation  developed  from  the  present 
heating  and  cooling  data  is: 

ft  = (1 . 25ln  (Ret)  - 0.235)  "2  (4) 

The  standard  deviation  of  the  correlated  ft  from  the  measured  value  was  estimated  to  be  1.9%  of  the 
measured  value. 

Figure  4 compares  the  present  measurements  of  ft  to  the  friction  factor  correlation  given  by  Filonenko 
(1954)  for  turbulent  flow  within  a smooth  tube: 

ft  = ( 1 . 58ln  (Ret)  - 3.28)'2  (5) 

The  Filonenko  (1954)  equation  agrees  within  approximately  10%  of  the  measured  tube-side  friction 
factor.  The  heating  and  cooling  friction  factors  nearly  coincide  for  Reynolds  numbers  above  20,000. 
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The  friction  factors  for  cooling  are  marginally  above  those  for  heating  for  Reynolds  numbers  below 
10,000;  the  difference  is  within  the  uncertainty  of  the  measurement  for  low  mass  flow  rates.  The 
agreement  between  the  present  measurements  and  the  Filonenko  (1954)  equation  indicates  that  the  method 
used  to  measure  the  friction  factor  is  reliable. 

The  measured  annulus-side  fanning  friction  factor  was  correlated  to  the  measured  annulus  Reynolds 
number.  The  resulting  annulus-side  fanning  friction  factor  correlation  for  the  spined  annulus  with  the 
spines  touching  the  outer  wall  is: 

fa  = (0 . 163ln  (Rea)  + 0.238)'2  (6) 

Figure  5 compares  the  measured  data  to  the  above  correlation.  The  estimated  standard  deviation  between 
the  measurements  and  the  correlation  was  2.4%  of  the  measurement.  No  measurable  difference  between 
heating  and  cooling  data  exists. 

The  friction  correlation  given  by  Hobson  (1954)  is  approximately  60%  below  the  data  presented  here. 
This  supports  Hobson’s  (1954)  conjecture  that  his  relationship  for  spines  in  contact  with  the  outer  wall 
would  not  be  satisfactory.  Apparently,  a significantly  lower  friction  factor  results  when  the  flow  bypasses 
the  fins.  The  measured  friction  factor  for  the  integral -spine-fin-annulus  is  15  times  the  laminar,  smooth 
tube  friction  provided  by  Davis  (1943).  This  is  expected  since  the  profile  drag  on  the  fins  contributes 
much  more  to  the  pressure  drop  than  does  the  skin  friction  on  the  tube  walls. 


HEAT  TRANSFER 


The  heat-transfer  coefficient  for  the  integral-spine-fin-annulus  (ha)  was  derived  from  the  measured  overall- 
conductance (UA)  of  the  heat  exchanger  using  a Wilson  plot  technique  similar  to  that  given  by  Briggs 
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and  Young  (1969).  The  methodology  of  the  present  Wilson  plot  differs  from  that  of  Briggs  and  Young 
(1969)  in  two  ways:  (1)  the  Prandtl  number  exponent  on  the  annulus  side  is  taken  as  unknown  and,  (2) 
an  additional  term  is  included  to  account  for  thermally  developing  flow  of  the  annulus  side. 


The  governing  equation  of  the  Wilson  plot  technique  is  the  overall  conductance  equation.  The  overall 
conductance  (UA)  can  be  represented  by  a sum  of  the  individual  thermal  resistances: 


UA  q hcAc  2n  Lkcu  r\haAs 


The  log-mean  temperature  difference  (ATlm)  is  calculated  from  the  measured  temperature  drop  of  each 
fluid  stream.  The  heat  load  (q)  of  the  heat  exchanger  is  calculated  from  the  measured  mass  flow  rate  and 
temperatures;  the  average  of  the  two  fluid  streams  is  used.  The  ht  is  the  heat-transfer  coefficient  of  the 
inner  smooth  tube  and  is  based  on  the  tube-side  area  At  (ttDjL).  The  next  to  the  last  term  is  the  thermal 
resistance  of  the  tube  wall  between  the  fluids.  The  only  parameters  in  Eq.  (7)  that  cannot  be  directly 
measured  or  calculated  are  the  tube-side  and  annulus-side  heat-transfer  coefficients,  ht  and  ha. 


The  basic  strategy  of  the  Wilson  plot  (1915)  is  to  substitute  functional  forms  with  unknown  constants  for 
the  tube-side  and  annulus-side  thermal  resistances  into  the  overall  conductance  equation  and  linearize  it 
so  that  the  unknown  constants  can  be  obtained  from  the  slope  and  the  intercept  of  a linear  regression. 
The  functional  form  of  the  tube-side  Nusselt  (Nut)  number  was  assumed  to  be: 


Nut 


h£, 

kc 


_1 

Ct  Reta  Pr  t3 


\ 0 . 14 


t 


(8) 


The  last  parameter,  the  ratio  of  the  bulk  and  wall  liquid  viscosities  to  the  0.14  power,  is  the  Sieder-Tate 
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(1936)  correction  to  account  for  radial  property  variations.  The  liquid  Reynolds  number  (Ret  = 
4mt/(/xt7rDi))  and  the  liquid  Prandtl  number  (Prt)  are  evaluated  at  the  bulk  temperature.  The  unknown 
constant  Ct  is  determined  via  the  Wilson  plot.  Conversely,  the  value  of  the  Reynolds  number  exponent 
(a)  is  specified  with  a certain  degree  of  uncertainty. 

The  dubiousness  of  the  tube-side  Reynolds  number  exponent  (a)  can  be  a serious  limitation  to  the  Wilson 
plot  method  (Shah,  1990).  Shah  (1985)  notes  that  the  exponent  (a)  may  vary  from  0.7  to  0.9  depending 
on  the  Reynolds  and  Prandtl  numbers.  To  remedy  this,  Webb  (1993)  suggested  that  the  exponent  (a)  be 
determined  from  a best  fit  of  Eq.  (8)  to  a modified  version  of  the  Petukhov-Popov  equation  (Gnielinski, 
1976)  — valid  in  the  range  2300  < Re  < 104  --  for  the  Pr  and  Re  range  of  our  data.  It  is  appropriate 
to  use  the  Petukhov-Popov  equation  to  determine  the  exponent  since  it  predicts  the  single-phase  Nusselt 
number  for  fully  developed  flow  within  a smooth  tube  to  an  uncertainty  of  ±6%  for  0.5<Pr<200.  A 
Reynolds  number  exponent  of  0.897  and  a Ct  of  0.0122  resulted  from  this  fit. 

The  methodology  of  the  Wilson  plot  dictates  that  one  must  assume  a functional  form  for  the  annulus-side 
Nusselt  number.  The  smooth  annulus  Nusselt  number  is  a constant  for  fully  developed  laminar  flow 
(Jakob  and  Rees,  1941).  However,  a substantial  proportion  of  the  total  test  section  length  was  in 
thermally  developing  flow  for  most  of  the  test  data.  Consequently,  a functional  form  is  necessary  to 
account  for  thermally  developing  flow.  Kays  and  Crawford  (1980)  provide  such  an  expression  for  a 
circular  tube  of  length  L and  inner  diameter  D j with  a constant  wall  temperature  (Tw): 


Re  Pr In 


(9) 


Nu. 


m 


It  is  speculated  that  Eq.  (9)  contains  the  primary  parameters  which  govern  thermally  developing  laminar 
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flow  in  a spine-fin-annulus.  The  product  of  Re,  Pr,  and  D;  is  also  found  in  the  expression  for  the  local 
Nusselt  number.  The  logarithmic  factor  (ln9)  results  from  the  integration  of  the  local  developing  Nusselt 
number  over  the  entry  length.  The  resulting  average  Nusselt  number  in  the  thermally  developing  region 
depends  on  the  entering  fluid  temperature  (Tj),  the  wall  temperature  (Tw),  and  the  bulk  temperature  of 
the  fluid  (Tm).  Average  values  of  Tw  and  Tm  were  used  since  the  wall  temperature  was  not  constant. 


Making  an  analogy  with  Eq.  (9),  the  following  form  was  chosen  for  the  annulus-side  Nusselt  number 
(Nua): 


hJDh 

Nua  = = CaReab  Pran 


{DhlnT"  " 
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T -T 


(10) 


The  above  form  produced  a Wilson  plot  with  the  best  fit  of  all  forms  attempted.  A best  fit  was  defined 
as  one  that  produced  the  smallest  percent  standard  deviation  ( % sY i and  %sYo)  from  the  Wilson  plot  fit 
with  a value  of  Ct  close  to  0.0122.  Table  2 shows  the  effect  of  the  inclusion  of  various  combinations  of 
the  Grashof  number  (Gr),  the  Sieder-Tate  (1936)  correction,  and  the  thermally  developing  correction 
(Dhpln0p/Lp)  proposed  here.  Only  the  Dhpln9p/Lp  term  was  included  in  the  annulus-side  correlation 
since  its  inclusion  without  the  other  corrections  produced  the  best  fit.  The  exclusion  of  the  Gr  and  the 
property  correction  suggests  that  natural  convection  and  large  radial  temperature  gradients  are  not  present. 
This  is  consistent  with  the  values  of  Gr/Re2  < < 1 (typically  0.015  for  our  tests),  and  Prandtl  numbers 
below  25.  The  Sieder-Tate  (1936)  correction  is  probably  not  necessary  since  for  the  given  Pr  and  Re 
ranges  the  fins  radially  mix  the  flow.  The  thermally  developing  term  was  a 4 to  20%  correction  of  Nua 
for  our  data. 

Substituting  the  tube-side  and  the  annulus-side  Nusselt  number  forms  into  the  overall  conductance 
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equation  and  rearranging  into  a form  that  is  convenient  for  the  linear  solution  of  Ct  and  Ca  gives: 
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Figure  6 shows  a plot  of  Eq.  (11)  where  the  ordinate  (Y j)  is  the  first  term,  and  the  abscissa  (Xj)  is  the 
second  term  divided  by  Ct_1.  The  constants  Ct  and  Ca  are  the  inverses  of  the  slope  and  intercept  of 
Figure  6,  respectively.  The  standard  deviation  of  the  measured  Yj  from  that  obtained  from  a linear  best 
fit  of  Eq.  (11)  is  2.8%. 

Equation  (11)  can  be  rearranged  into  a convenient  form  for  the  solution  of  the  Reynolds  number  exponent 
(b)  as: 
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b lnRea  + In (Ca) 


The  first  term  (Y2)  of  Eq.  (12)  is  plotted  versus  ln(Rea)  which  is  X2  in  Figure  7.  The  annulus  Reynolds 
number  exponent  was  obtained  from  the  slope  of  a linear  best  fit  of  Eq.  (12).  The  constant  Ca  was 
obtained  from  the  intercept  of  this  fit.  The  Ca  from  Eq.  (12)  must  agree  with  that  obtained  from  Eq. 
(11).  The  standard  deviation  of  the  measured  Y2  from  that  obtained  from  a linear  best  fit  of  Eq.  (12) 
is  3.8%. 
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The  additional  unknown  exponents  n and  p were  solved  iteratively.  First,  values  for  n and  p were 
assumed.  Then  Ct  and  Ca  were  obtained  from  a linear  regression  of  Eq.  (11).  The  resulting  Ct  was 
substituted  into  Eq.  (12)  to  solve  for  an  updated  annulus-side  Reynolds  number  exponent  (b)  and  Ca. 
Iterations  on  Eqs.  (11)  and  (12)  were  performed  until  successive  values  for  Q and  Ca  converged.  An 
updated  value  for  the  annulus-side  Prandtl  number  exponent  (n)  and  the  exponent  p were  obtained  from 
a regression  on  the  annulus-side  Nusselt  numbers.  The  annulus  Nusselt  numbers  were  solved  by 
substituting  the  tube-side  heat-transfer  coefficient  obtained  from  the  correlation  developed  from  the 
Wilson  plot  into  the  overall  conductance  equation.  The  procedure  was  repeated  until  all  constants 
converged  to  a single  value. 

The  Ct  from  the  Wilson  plot  was  0.0122.  The  tube-side  Nusselt  number  obtained  from  the  Wilson  Plot 
agreed  with  the  Petukhov-Popov  equation,  on  average,  within  3.3%. 

The  surface  efficiency  (77)  was  approximated  by  assuming  that  the  fin  heat-transfer  coefficient  was 
constant  and  the  fin  cross  section  was  uniformly  rectangular.  The  fin  cross  sectional  area  and  perimeter 
were  estimated  from  length  weighted  average  values  of  the  tail  and  head  of  the  fin  of  the  spine  shown 
in  Figure  2.  The  average  77  was  0.74.  The  fin  efficiency  calculation  procedure  developed  by  Carranza 
(1993)  for  pyramid  shaped  spines  was  not  used  since  it  was  felt  that  the  fins  in  this  study  more  closely 
resembled  rectangles  rather  than  pyramids. 


Nu„  Correlations 

cl 

The  following  empirical  Nusselt  number  correlation  for  the  integral-spine-fin  annulus  was  obtained  from 
the  Wilson  plot: 
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Equation  13  is  recommended  for  Prandtl  numbers  from  4 to  22,  and  Reynolds  numbers  from  100  to  1400 
where  Gr/Re2  < < 1.  The  uncertainty  of  Eq.  (13)  was  estimated  from  a propagation  of  errors  to  be, 
on  average,  17%  of  Nua. 

The  solution  of  Eq.  (13)  requires  an  iteration  for  the  thermal  developing  correction.  If  an  explicit 
solution  is  preferred,  a version  of  Eq.  (13)  without  the  thermal  developing  correction  is: 

Nua  = 0 . 142Rea0,508  Pra°'416  (14) 

The  uncertainty  of  Eq.  (14)  is  the  same  as  that  of  Eq.  (13).  However,  a small  Nua  uncertainty  does  not 
verify  the  form  of  the  Nusselt  number  equation;  it  originates  from  small  uncertainties  in  the  calculation 
of  Yb  X1?  Y2,  and  X2.  A good  Wilson  plot  fit  --  small  sYi  and  sY2  — suggests  that  the  chosen  form  is 
correct  for  the  correlated  data  set.  If  the  genuine  form  was  known,  the  correlation  would  predict  all 
ranges  of  data  to  within  its  calculated  uncertainty.  Moreover,  a correlation  could  have  a small  uncertainty 
and  be  unable  to  predict  outside  the  data  source  range  because  it  has  an  incorrect  form.  Therefore,  the 
authors  prefer  Eq.  (13)  over  Eq.  (14)  since  it  is  expected  that  errors  may  result  using  Eq.  (14)  for  heat 
exchangers  that  do  not  experience  thermally  developed  flow  similar  to  that  experienced  by  our  test 
section. 

A comparison  of  the  exponents  in  Eqs.  (13)  and  (14)  to  exponents  of  other  empirical  correlations  in  the 
literature  somewhat  authenticates  Eqs.  (13)  and  (14).  First,  the  Reynolds  number  exponent  is  between 
that  obtained  by  Chen  et  al.  (1946)  for  laminar  flow  (0.45)  and  that  obtained  by  Davis  (1943)  for 
turbulent  flow  (0.8)  within  a smooth  annulus.  Even  though  the  present  annulus  Reynolds  numbers  are 
laminar,  the  staggered  fins  presumably  produce  chaotic  mixing  which  is  more  effective  than  laminar  flow. 
This  mixing  is  not  quite  as  efficient  as  the  random  mixing  of  turbulent  flow.  Consequently,  the  Reynolds 
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number  exponent  reflects  an  intermediate  value  between  turbulent  and  laminar  flow.  Second,  Chen  et 
al.  (1946)  used  (Dh/L)0  4 to  account  for  the  thermally  developing  boundary  layer.  The  power  to  which 
the  Dh/L  term  is  raised  in  the  Chen  et  al.  (1946)  correlation  is  close  to  the  power  used  for  the  same  term 
in  Eq.  (13). 


The  influence  of  the  spines  on  the  enhancement  of  the  heat  transfer  can  be  investigated  with  the  Chen  et 
al.  (1946)  correlation  which  was  developed  for  laminar  flow  within  an  annulus  without  fins: 


Nu_  = 1.02  Ref ' 45  Pr  =3| 
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For  example,  if  the  spine-fin  hydraulic  diameter  is  used  for  Dh  along  with  the  pre-finning  pipe  diameter 
Dp  for  Dr,  Eq.  (15)  underpredicts  the  integral-spine-fin  Nusselt  number  by  an  average  of  26%.  This 
suggests  that  the  heat  transfer  enhancement  is  due  not  only  to  the  additional  surface  area,  but  also  to  the 
flow  mixing  that  the  fins  provide.  Consequently,  good  heat  transfer  practice  dictates  that  the  fins  touch 
the  outer  wall.  Otherwise,  as  shown  by  Knudsen  and  Katz  (1950),  laminar  flow  remains  in  unmixed 
ribbon-like  stream  lines  in  the  annular  space  between  transverse  fins  and  the  outer  wall. 

The  magnitude  of  the  fins  enhancement  can  be  better  illustrated  with  a comparison  of  the  fin  and  root  heat 
transfer  coefficients.  The  fin  efficiency  weighted  heat-transfer  coefficient  for  the  fin  (rjfhf)  normalized 
by  that  for  the  root  surface  (hr)  is  plotted  versus  the  Prandtl  number  in  Figure  8.  The  hr  was  obtained 
from  the  Chen  et  al.  (1946)  correlation  using  the  root  diameter.  The  hf  was  obtained  from  substitution 
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of  hr  and  the  measured  spine  fin  heat-transfer  coefficient  into  a sum  of  the  individual  heat-transfer 
coefficients  weighted  by  surface  area: 


*f= 


(16) 


The  figure  shows  that  the  fin  has  anywhere  from  10%  to  50%  higher  heat  transfer  that  the  root  surface. 
Presumably,  the  fin  enhancement  results  from  high  heat  transfer  at  the  leading  edge  of  the  fins 
(performing  like  short  plates)  and  the  flow  mixing  produced  by  the  staggered  fins.  The  mixing 
enhancement  becomes  more  pronounced  at  higher  Pr  since  the  fins  are  more  effective  than  the  root 
surface  at  converting  momentum  to  thermal  energy. 


Extrapolation  of  Eq.  (13)  for  spined  annuli  with  different  hydraulic  diameters  was  done  by  including  the 
same  parameter  ((D0/Dh)°  8)  that  Chen  et  al.  (1946)  used  to  account  for  the  effect  of  different  smooth 
annuli  geometries.  The  resulting  equation  is: 
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Equation  (17)  should  not  be  used  if  the  heat  transfer  is  governed  by  natural  or  mixed  convection. 
Moreover,  the  extension  of  Eq.  (13)  to  other  hydraulic  diameters  is  only  an  intuitive  guess;  data  for 
annulus  spine-fin  of  various  hydraulic  diameters  are  required  to  verify  or  disprove  it. 


CONCLUSIONS 


Empirical  correlations  for  the  laminar,  single-phase  heat  transfer  and  friction  characteristics  of  an  integral- 
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spine-fin  pipe  within  an  annulus  have  been  presented.  The  heat-transfer  coefficient  was  determined  with 
a modified  version  of  the  Wilson  Plot  method.  An  additional  term  was  included  in  the  empirical 
correlation  for  the  Nusselt  number  to  account  for  the  development  of  the  thermal  boundary  layer.  The 
exponential  dependence  of  this  term  and  the  annulus-side  Reynolds  and  Prandtl  numbers  were  all 
determined  via  the  modified  Wilson  plot.  The  single-phase,  smooth  tube-side  correlation  obtained  from 
the  Wilson  plot  agreed  to  within  3%  of  the  Petukhov-Popov  equation.  The  spines  enhance  the  heat 
transfer  through  additional  surface  area,  fluid  mixing  and  performing  similar  to  short  plates. 
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Table  1 Average  Estimated  Percent  Measurement  Uncertainties 
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Table  2 Best-Fit  for  Annulus-Side  Nusselt  Number 
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Figure  1 Spine-fin  surface  and  heat  exchanger  geometry 
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Figure  2 Schematic  of  spine 
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Figure  3 Schematic  of  test  rig 
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Figure  4 Measured  smooth  tube  friction  factor  (ft) 
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Figure  5 Measured  integral-spine-fin  annulus  friction  factor  (fa) 
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Figure  6 Wilson  plot  1 
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Figure  7 Wilson  plot  2 
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Figure  8 Comparison  of  the  relative  magnitudes  of  the  fin  heat-transfer  coefficient  to  that 
for  the  root  surface 
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